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^~~l Abstract. We study analytically a variant of the one-dimensional majority-vote 

^_^ model in which the individual retains its opinion in case there is a tie among the 

,^ neighbors' opinions. The individuals are fixed in the sites of a ring of size L and 

can interact with their nearest neighbors only. The interesting feature of this model 
is that it exhibits an infinity of spatially heterogeneous absorbing configurations for 
I L — > oo whose statistical properties we probe analytically using a mean-field framework 

jrt based on the decomposition of the L-site joint probability distribution into the n- 

contiguous-site joint distributions, the so-called n-site approximation. To describe 
the broken-ergodicity steady state of the model we solve analytically the mean-field 
dynamic equations for arbitrary time t in the cases n = 3 and 4. The asymptotic 
Ch limit t — > cx) reveals the mapping between the statistical properties of the random 

i-rt initial configurations and those of the final absorbing configurations. For the pair 

^H approximation {n = 2) we derive that mapping using a trick that avoids solving the full 

dynamics. Most remarkably, we find that the predictions of the 4-site approximation 

reduce to those of the 3-site in the case of expectations involving three contiguous sites. 

.. In addition, those expectations fit the simulation data perfectly and so we conjecture 

^ that they are in fact the exact expectations for the one-dimensional majority-vote 

CN| model. 

\o 
o 

^ PACS numbers: 89.65.-s, 89.75.Fb, 87.23. Ge, 05.50.-Fq 

O 

(N 

T-H 1. Introduction 

> 

K^ A desirable property of a model for social behavior, or for complex systems in general, 

;-( is the presence of a nontrivial steady state characterized by infinitely many equilibrium 

points in the thermodynamic limit. This was the main appeal of the mean-field spin- 
glass models used widely since the 1980s to study associative memory pTj, prebiotic 
evolution [2], ecosystem organization [3], and social systems |1] to name just a few of 
the areas impacted by the spin-glass approach to model complex systems [5]. 

Models of social dynamics, however, are typically defined through the specification 
of the dynamic rules that govern the interactions between agents [B] and so they are not 
amenable to analysis using tools borrowed from the equilibrium statistical mechanics 
of disordered systems. Nevertheless, the display of a steady state characterized by 
a multitude of locally stable and spatially inhomogeneous configurations remains a 
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celebrated feature of this class of models, whose paradigm is Axelrod's model [7], 
since they can explain the diversity of cultures or opinions observed in human societies. 
Axelrod's model is attractive from the statistical physics perspective because it exhibits 
a nonequilibrium phase transition which separates the spatially homogeneous (mono- 
cultural) from the heterogeneous (multicultural) regimes [5| [^ [TU]. 

More recently, a long-familiar model of lattice statistical physics - the majority- 
vote model [m |T2] - was revisited in the context of social dynamics models [131 [T3] . In 
fact, the majority-vote model is a lattice version of the classic frequency bias mechanism 
for opinion change [I9j, which assumes that the number of people holding an opinion 
is the key factor for an agent to adopt that opinion, i.e., people have a tendency to 
espouse opinions that are more common in their social environment. The variant of 
the majority-vote model considered in those studies includes the state of the target 
site (the voter) in the reckon of the majority (hence we refer to the model as extended 
majority- vote model), which happens to be the variant originally proposed in the physics 
literature [HI [I2]. This fact is the sole responsible for the existence of an infinity 
of heterogeneous absorbing configurations whose statistical properties were thoroughly 
studied via simulations in the case of a two-dimensional lattice [13] . Moreover, the non- 
linearity of the transition probabilities resulting from the inclusion of the voter opinion 
in the majority reckoning makes the model not exactly solvable, in contrast to the voter 
model for which the transition probabilities are linear [15] . 

Many interesting variants of the one- dimensional majority-vote model were 
considered in the literature. For instance, some variants separate the individuals in 
groups of fixed sizes and apply the majority-vote rule to update the opinion of the 
entire group simultaneously [THj. Others differentiate the groups a priori by introducing 
a group-specific bias used to determine the group opinion in case of ties [17J. Another 
variant of interest is the non-conservative voter model for which the probability that the 
voter changes its opinion depends non-linearly on the fraction of disagreeing neighbors 
[18j . In particular, this variant reduces to the model we study in this paper in the 
case the voter changes opinion solely when confronted by a unanimity of opposite- 
opinion neighbors. All these variants of the majority-vote rule model were studied 
via simulations or within the single-site mean-field framework, except for the non- 
conservative variant which was examined within the pair approximation as well [T8] . 
Here we show that the single-site and the pair approximations yield incorrect predictions 
for all statistical measures of the steady states and argue that the 3-site and 4-site 
approximations yield the exact results for measures involving up to three contiguous 
sites of the chain. 

In this contribution we study analytically the one- dimensional version of the 
extended majority- vote model, which is described in Sect. [2| Our goal was to understand 
how the multiple-cluster steady state of the model could be described within the mean- 
field approach (Sect. [s]). We find that the signature of the ergodicity breaking is the 
appearance of an infinity of attractive fixed points in the mean-field equations for the 
n-site approximation with n > 2. The characterization of the mean-field steady state 
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requires the complete analytical solution of the dynamics in order to obtain the mapping 
between the statistical properties of the random initial configurations and those of the 
final absorbing configurations, except for the 2-site or pair-approximation for which we 



find a simple shortcut to that mapping, as described in Sect. |3.2[ The full solution 
of the dynamics was obtained for the 3 and 4-site approximations in Sects. 3^ and 
3.4[ respectively. We find that these two approximation schemes yield the very same 
expressions for the steady-state expectations involving three contiguous sites [see Eqs. 



(24), (34) and (35)] and so we conjecture that those expressions are exact. A perfect 
fitting of the simulation data by those predictions adds further support to this claim. 
In addition we find that the steady-state expectations involving four contiguous sites 
calculated within the 4-site approximation fit the simulation data perfectly. However, 
this approximation fails to describe higher order expectations. 

2. Model 

The agents are fixed at the sites of a ring of length L and can interact with their 
nearest neighbors only. The initial configuration is chosen randomly with the opinion 
of each agent being specified by a random digit 1 or with probability po and 1 — po, 
respectively. At each time we pick a target agent at random and then verify which is 
the more frequent opinion (1 or 0) among its extended neighborhood, which includes 
the target agent itself. The opinion of the target agent is then changed to match the 
corresponding majority value. We note that there are no ties in the calculation of the 
preponderant opinion since the extended neighborhood of any agent comprises exactly 
3 sites. As a result, the update rule of the model is deterministic; stochasticity enters 
the model dynamics through the choice of the target site and in the selection of the 
initial configuration. The update procedure is repeated until the system is frozen in an 
absorbing configuration. 

Although the majority-vote rule or, more generally, the frequency bias mechanism 
for cultural change [19] is a homogenizing assumption by which the agents become 
more similar to each other, the two-dimensional version of the above-described model 
does exhibit global polarization, i.e., a non-trivial stable multicultural regime in the 
thermodynamic limit [131 HZl IH]. This regime should exist in the one-dimensional 
version as well, since any sequence of two or more contiguous I's (or O's) is stable 
under the update rule. It should be noted that for the more popular variant of the 
majority- vote model, in which the state of the target site is not included in the majority 
reckoning, and ties are decided by choosing the opinion of the target agent at random 
with probability 1/2, the only absorbing states in the thermodynamic limit are the 
two homogeneous configurations [201 121] . As mentioned before, the inclusion of the 
target site in the calculation of the majority is actually the original definition of the 
majority- vote model as introduced in Refs. [II1[I2|- Figure [I] illustrates an absorbing 
configuration of the extended majority- vote model together with a random configuration 
with the same density of I's. The larger number of clusters (domains) observed in the 
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Figure 1. (Color online) Disposition of the Ui = 1 variables in a ring with L = 500 
sites. The inner circle (black) is an absorbing configuration of the extended majority- 
vote model with p — 0.49. There are 92 clusters and the largest cluster comprises 19 
sites. The outer circle (blue) shows a random configuration with the same density of 
I's. The total number of clusters is 255 and the largest one comprises 9 sites. 



random configuration is due to the possibility of isolated sites, which are unstable under 
the majority- vote rule. 

As usual, we represent the state of the agent at site i of the ring by the binary 
variable ai = 0, 1 and so the configuration of the entire ring is denoted by a = 
(cTi, (72, . . . , ctl). The master equation that governs the time evolution of the probability 
distribution P (a, t) is given by 

|P (a, t) = J2 [W^ (^0 P i^'^t) - W. (a) P {a, t)] (1) 

i 

where a* = (cxi, . . . , 1 — cXi, . . . , o"l) and Wi (a) is the transition rate between 
configurations a and a* [201 |2T]. ^or the one- dimensional extended majority- vote model 
we have 

Wi (cr) = (Tj (1 - CTj-l - CTj+l) + (Ti^lCTi+l (2) 

for i = 1, . . . ,L. The boundary conditions are such that (Tq = a^ and ctl+i = o"i. To 
implement the n-site approximation up to n = 4 we need to evaluate the following 
expectations 

|(a.) = ((l-2a,)W^,(a)), (3) 

|(a.a,) = 2(a,(l-2a,)iy,(a)), (4) 
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— (aiajak) = 3 {ajak (1 - 2(7^) Wi (a)) , (5) 



{a-iajakai) = 4 (o-^crfca/ (1 - 2cri) Wi (a)) (6) 



d 
dt 

where all indexes are assumed distinct and we have introduced the notation ((...)) — 
So- (• • •) -P ('^5^)- The n-site approximation is based on the calculation of this average 
by replacing the full joint distribution probability P {a, t) by a decomposed form that 



depends on the order n of the approximation [see Eqs. (|8]), (11), (18) and (42)]. Of 
course, in the derivation of Eqs. ([3])-(|6]), which generalize trivially to an arbitrary number 
of sites, we have assumed translational invariance, i.e., all sites are assumed equivalent. 

3. Mean-field Analysis 

In this section we study the one-dimensional extended majority-vote model using the 
well-known mean-field n-site approximation (see [221 ISSl [211 [25] ) . The basic idea behind 
the n-site approximation is to rewrite the distribution P {a, t) in terms of elementary 
joint probabilities of n contiguous sites only and then deriving a system of self-consistent 
equations for these probabilities. This key idea is expressed mathematically using the 
following equation which summarizes the approximation scheme, 

Pl\L-l (CTJ I (Ti, . . . , (Tl) = Pl\2n-2 {(^i \ CTj-n+l, • • • , (^i+n-l) (7) 

where, of course, 0", does not appear in the arguments at the right of the | delimiter in 
these conditionals. This procedure will be illustrated in the next subsections for n = 1 to 
n = 4. It is interesting to note that, except for the single-site approximation n = 1, the 
states of any two sites are statistically dependent variables regardless of their positions 
in the ring. 

3.1. The single-site approximation 

This is the simplest mean-field scheme which assumes that the state of the agents at 
different sites are independent random variables so that 

P {a, t) = pi (cTi, t) pi ((72, t)---Pi (o-L, t) (8) 

and so it is only necessary to calculate the one-site distribution pi {ai, t) to describe the 
dynamics completely. This can be done by noting that p = {cyi)t = Pi (1,^) and using 
Eq. (pi) to derive a self-consistent equation for p. The final result is simply 

p = p(-2p2 + 3p-l) (9) 

with the notation x = dx/dt. We note that p contains the same information as the 
single-site probability distribution pi (ai) since pi {ci = 1) = p and pi (cXi = 0) = 1 — p. 
A straightforward stability analysis shows that there are three fixed points, 

Pl = 0, P2 = l/2, P3 = l, (10) 
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with only pi and ps being stable. This means that only the homogeneous configurations 
are stable and so the single-site approximation completely fails to describe the steady 
state of the extended majority-vote model. 

3.2. The pair approximation 

Using Eq. ([T]) to write the full probability distribution in terms of the joint probability 
of two sites and omitting the time dependence we find 

^(cr) = J — ^ J — ^ J ^ J — ^ (11) 

Pi (o"i) Pi {(^2) ■■■Pi {aL-i)Pi [aL) 

where pi (ctj) = ^o-,P2 (o"j, ctj)- To avoid overburden the notation we use the same 

notation for pi as done in the single-site approximation but the pi which appears in Eq. 



(11) is numerically distinct from that calculated in the previous section. This simplifying 
convention for the notation of probabilities will be used in the next sections as well. 

Within this framework, it is only necessary to calculate p2 (o"j5cri+i) to describe 
the dynamics of the model completely. This amounts to 4 variables, namely, 
P2 (1, 1) ,P2 (1, 0) ,p2 (0, 1) and p2 (0, 0), but use of the normalization condition and of 
the parity symmetry [p2 (1, 0) = p2 (0, 1)] allows us to reduce the number of independent 
variables to only 2. The first variable we pick is = (ctjCTj+i) = p2 (1, 1) which is given 
by Eq. ^ with j = i-\-l. Next, noting that p2 (1, 0) = pi (1) — p2 (1, 1) = p — <P we pick 
p, given by Eq. ([3]), as the second independent variable. Carrying out the averages in 
the right-hand sides of Eqs. ^ and ^ using the decomposition (11) yields 

\2 



(p-0)^(2p-l 
P = 

and 



'^- 2;(i-;) <^'' 



(p- 



The steady-state condition = p = as well as the numerical integration of these 
equations yields p = cj) for t — )■ oo, with p determined by the value of the initial 
condition p (t = 0) = po and (t = 0) = Pq. We note that this result implies that 
P2 (1, 0) = p2 (0, 1) = 0, meaning that the number of interfaces between clusters, i.e., of 
contiguous sites in different states at the steady state, is not extensive. This prediction is 
not correct as indicated by the higher-order approximations and by the simulation data. 
Despite this incorrect prediction, the pair approximation explains the most remarkable 
aspect of the extended majority-vote model, namely, the ergodicity breaking reflected 
by the inflnity of distinct absorbing conflgurations. 

The imposition of the steady-state condition is not sufficient to determine the 
equilibrium solution p = (f) because there is a continuum of fixed points characterized 
by the function p (po). In order to obtain this function or map, we revert to the original 



variable x = P2{1,0) = p — 0, and rewrite Eqs. (12) and (13) as 
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X 



X 



2p(l-p)- 

from which we can immediately obtain the integral equation 

dp' 



x{t) 



dx' 



Pit) 



1' 



(15) 



(16) 



As the stationary regime is obtained in the limit t — > oo, we define p (po) = p{t ^ oo). 
In addition, using xq = po (1 — Po) and a; (t — t- oo) = (steady-state condition) we find 

1 



p(Po) 



[l + (2po-l)e 



2po{l-po) 



]. 



(17) 



This equation is identical to that derived for the non-conservative voter model in the 
case the voter changes opinion only when confronted by a unanimity of opposite-opinion 
neighbors [IS]. Figure |2] shows this steady-state solution together with the results of the 
simulations for a ring with L = lO'' sites. Despite the incorrect prediction (0 = p), Eq. 



(17) yields a remarkably good quantitative agreement with the density p obtained from 



the simulations. However, as we will show next the (supposedly) exact expression for p 



is much simpler than Eq. (17) 




Figure 2. (Color online) The fraction of sites in state 1 at equilibrium, p (filled 
circles), and the probability that two neighbors are in state 1 at equilibrium, (/> (filled 
triangles) as functions of the initial fraction of sites in state 1, po- The solid line is 
the result of the pair approximation for which p = (j) given by Eq. (17). The initial 



condition is p = po a-nd (p ^ p^. The symbols show the results of the simulations for a 
ring of size L = 10^ and 10^ independent samples. 



3.3. The 3-site approximation 

In this scheme the decomposition of P (a) is 

P3 (C^l, 0-2, 0-3) P3 (0-2, CTS, (T4)... Ps {aL, CTi, CTs) 



Pia) 



P2 (CTI, (T2) P2 (cr2, 0-3) .. . P2 {ai, CTl) 



(18) 
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where p2 (o"i,crj+i) = J2a Ps (c^j, ctj+i, aj+2)- The goal here is to calculate the 9 
probability values ps (ctj, ctj+i, cri+2) with dk = 0, 1. As before, use of the normalization 
condition and of the parity symmetry give us 6 variables to be determined using 
appropriate linear combinations of Eqs. (|3])-([6]). We choose the following variables 

a^o =P3 (0,0,0), xi =1)3(1,0,0), X2 =P3 (1,1,0), 

a^ic = P3 (0, 1, 0) , X2C = Pa (1,0,1), X3 = ps (1, 1, 1) . 

which are given by the expectations Xi = {ai) — (ctjCTj+i) — {o-iO-i+2) + (o"'jO"'j+iO"i+2), 
X2 = ((XiCTj+i) - (o-iCri+i(Ti+2), 3:3 = (o-iai+iorj+2), and so on. 

3.3.1. Mean-field equations. Evaluating the averages in Eqs. ([3])-(|6| using the 



decomposition (18) yields 



1 Xic / X .1 X2C , X 

a;i = ■ \X2c-X\)-, ^■i=o . [X1C-X2), 

6 Xic + X2 3 X2C + Xi 

1 Xic 1 X2C 

xic = -7; \ (3x2c + a;i), X2c = -7; \ (3xic + 0:2) , (20) 

3 X2C + Xi 3 Xic + X2 

1 Xic 1 X2C 

xq = ^ [xic + 2x1 + X2) , X3 = ■ {X2c + xi + 2x2) . 

3 Xic + X2 3 X2C + Xi 

The steady state is given by xic = X2c = 0, i.e., ps (0, 1, 0) = ps (1, 0, 1) = which, in 
contrast to the situation we found in the analysis of the pair approximation, reflects the 
physical requirement that absorbing configurations cannot exhibit isolated sites. As we 
are still left with 4 undetermined variables after applying the steady-state condition, we 
need an alternative method to characterize the steady state. Somewhat surprisingly, in 
this case we will be able to solve the dynamics analytically, a feat that seems unfeasible 
in the case of the pair approximation. 



In fact, what makes the system of nonlinear coupled equations (20) solvable is the 
observation that y = xic + X2 = X2c + xi = p2 (1, 0), so the denominators in the r.h.s. 
of all those equations are identical. In addition, we note that xq does not affect the 
other 5 variables so we can first solve for them and then return to the equation for xo 



to complete the solution of the system ( 20 ) . 



Introducing the auxiliary variables zi = xic + X2c, ^2 = xic — X2c, and recalling 



that p = Xi -\- X2 + X2C + X3 we reduce Eqs. (20) to 



^1 . Zi Zi — Z2 . Z2 . Z2 , . 

where we have omitted the equation for xq. The last two equations can be immediately 
solved and yield 

2;2(t) = po(l-po)(l-2po)e-5*, (22) 

p it) = p^ (3 - 2po) + po (1 - po) (1 - 2po) e-5*, (23) 
Hence in the asymptotic limit t — ?■ 00 we obtain 

P(Po)=Po(3-2po) (24) 
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and Z2 = 0, as expected, since both xic and X2c vanish at the steady state. For po ^ 



or po ~^ 1 the pair approximation estimate of p given by Eq. (17) reduces to Eq. (24) in 
a first order approximation in po or 1 — po- Equation (24) describes the simulation data 



perfectly as illustrated in Fig. |3]and, as already mentioned, we believe it gives the exact 
value for the steady-state density of Is of the one- dimensional majority-vote model. 




Po 

Figure 3. (Color online) The solid lines are the analytical results for the steady- 
state measures obtained with the 4-sitc approximation while the symbols represent the 
results of the simulations for a ring of size L = lO"* and 10^ independent samples. The 
convention is (from top to bottom) p{po) (circles), (j){po) (triangles), ip (po) (squares) 
and wi (po) (upside down triangles). The upper three curves are identical for the 3-site 
approximation. 



The explicit calculation of the remaining two unknowns y and zi using Eqs. (21) 



is not too involved and their knowledge will allow us to evaluate other quantities of 
interest, such as and other high-order correlations. We begin by introducing the 
auxiliary variables ui = zi/y and U2 = Z2/y which satisfy the equations 



iOi 



-Ui 



-bJ. 



2) 



1 1 

UJ2 = --UJ2+ ::UJlUJ2- 



(25) 



3 3 

Next, we define a = wj — u;| which is given by d = — 2a/3 and so 

2t 

a (t) = aoe~^, (26) 

with cto = 4po (1 — Po)- At this point we can readily write an explicit equation for zi in 
terms of y, 



z^{t) = e ^^z^ + aoy^{t), 



(27) 



where we used the fact that zi >0 and that Z2 is given by Eq. (22). Now, inserting this 



expression in the equation for y [see Eq. (21)] results in an easily solvable integral, 



yo 



dy' 



log 



'nin 



y. 



yo+ V^o + vl 



(28) 
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with ml = Zq/oq and yo = po{l — po). Finally, carrying out the integration yields 
y{t) = ^Ivo + yml + yljexp 

(29) 



rrin 



exp 



2^/0 + V^o + Vl 
In the limit t — > oo this equation reduces to 



oo ( 1 — e~3 



y{po) =Po(l -Po)cosh 2a/po(1 -po) 
1 



Po (1 - Po) sinh 



2Vpo(l-po) 



(30) 



which exhibits the symmetry y (po) = y{l — po). 



To conclude the solution of the system of equations (20) we need now to determine 



xq. The easiest way to do that is to rewrite the equation for xq in (20) as 



Xq = -2y + - (ii - Sia) 
which can be immediately integrated to yield 

xo (t) = xo (0) - 2 [y (t) - z/o] + ^ h {t) - zi 
In the asymptotic limit xq (t — !■ oo) = -i/^-i we find 

V'-i(po) = (l + 2po)(l-po)'-2y(po). 



(31) 



- - [Z2 (t) - Z2 (0)] . (32) 



(33) 



3.3.2. Simple steady-sate expectations. At this stage we should be able to express any 
quantity characterizing the ring configuration at time t in terms of the time-dependent 
variables p, Z2, y and Zi. However, we will focus here only in the steady-state regime 
(t — 7- oo) for which only p and y contribute since Zi = Z2 = 0. 

The most interesting expectations are those whose time evolution are defined by 
Eqs. ([3])-(|6]) in the case the indexes are associated to contiguous sites. We begin with 
= ((TiO"j+i) = X2 + X3, and recall that p = (ctj) = xi + X2 + X2c + ^3 and y = xi -\- X2c 
so that (p = p — y. Then at the steady state we find 

0(Po)=Po(3-2po)-y(po). (34) 

Next we note that ip = (o"jO"i+io'j+2) = x^ = p — 2y + {zi + Z2) /2 and so 

V^(Po)=Po(3-2po)-2y(po). (35) 

These two steady-state expectations, which are shown in Fig. [3| describe the simulation 
data perfectly. Expectations involving more than three contiguous sites must be 
decomposed so as to be described by the 3-site approximation. Of particular interest 
is the 4-site expectation wi = (o-icri+io-i+20'j+3) = p4 (1,1, 1,1), which in the 3-site 
approximation scheme becomes Wi = pi (1, 1, 1) /p2 (1, 1) so that Wi (po) = ip^/(p at the 
steady state. In the scale of fig. [3] this result is indistinguishable from the simulation 
data or from their counterpart calculated with the 4-site approximation (see fig. 111). 
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11 



For completeness, let us calculate 0_i = p2 (0, 0) = xi + xq at the steady state. 
Since Xq is given by Eq. (33) and Xi = y — {zi — Z2) /2 we find 

0_i (po) = (1 + 2po) (1 - Pof - y (Po) . (36) 

The fact that the dynamics is invariant to the interchange of Is and Os provided that we 
change po to 1 — po is expressed by the easily verifiable identities p (1 — Po) = 1 — P (Po)) 
(po) = 0-1 (1 — Po), and -0 (po) = "^-i (1 — Po)- Of course, this symmetry holds for all 
orders n of the n-site approximation scheme and we will resort to it to abbreviate the 



calculations of the 4-site approximation in Sect. 3.4 



3.3.3. Probability of clusters of length m. A more informative quantity is the 
probability of finding a cluster of m > 1 sites in an absorbing configuration. There 
are only two possibilities for such a cluster: (a) a site in state Uj = followed by 
m sites in states cxj+i = o"j+2 • • • cn+m = 1 which are then followed by another site 
in state cTj+m+i = and (b) a site in state ctj = 1 followed by m sites in states 
cTj+i = (Ji+2 • • • CTj+m = which are then followed by another site in state (jj+m+i = 1- 
The probability of these configurations happening in an absorbing configuration can be 



easily derived using the decomposition ( 18 ) and yields 

-2 



Pf (m) 



P^ (0,1,1) 
P2(l,l) 



P3 (1, 1, 1) 



+ 



pj (1,0,0) 

P2(0,0) 



P3 (0,0,0) 

L P2(o,o) J 



.(37) 



P2(l,l) J 

To rewrite this expression in terms of more elementary steady-state quantities we recall 
that P3 (1, 0, 0) = p2 (0, 0) - p3 (0, 0, 0) and ^3 (0, 1, 1) = P2 (1, 1) - Ps (1, 1, 1) so that 



>(3) 



P^i'iPo,m) 



-ip) fip 



m-2 



+ 



^'-i-V^-i) ('ip 



m-2 



(38) 



The 3-site approximation estimate for the probability of finding clusters of length 
m given by this equation is presented in Figs. |5] and [6] together with the results of 
the simulations and the estimate of the 4-site approximation. We will postpone the 
discussion of the physical implications of the results presented in these figures to Sect. 

SI 



3.3.4. Two-site correlations. Knowledge of the two-site correlations defined by 

corr ((Ti, ai+j) = {(Xiai+j) - {ai){ai+j) (39) 

is very useful to determine the validity of the approximations. Since all sites are 
equivalent we have {(Xi) = (ctj+j) = p(po) regardless of the order n of the approximation. 
Some two-site expectations follow straightforwardly from the previous results, namely, 
{o'i<Ji)(3) = p, (crjcri+i)(3) = 4> and {<Ji<Ji+2) (3) = i'- The first nontrivial two-site 
expectation is 



((Tj(Ti+3)(3) = Z^ P((Tj = l,Crj+i,(Tj+2,Crj+3 = 1) 

_r , (0-1-^-1)' 



(40) 
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where we have decomposed the 4-site probabihty in terms of the elementary 3-site 
probabihties. Note that the sum has two non-vanishing terms only, since P (1, 0, 1, 1) = 
P(l, 1,0, 1) = 0. Applying the very same procedure to calculate (o"iO"j+4) yields 

In this case only 4 terms give nonzero contributions to the sum over the middle sites. 



Equations (40) and (41) clarify a fact that is often unappreciated, namely, regardless of 
their position in the ring, the sites are always treated as statistically dependent variables 
within the n-site approximation scheme for n > 1. 

3.4- The 4-site approximation 

In the 4-site approximation framework the decomposition of P (a) is given by the 
prescription 

„. ^ P4(ai,a2, (73,^4)^4(^2, as, (14, as)... P4(cxL,cxi,fT2,a3) 

^ (^) = 7 N ? N 7 N (42) 

P3 (ai, a2, as) ps (a2, as, a4) . . . ps (ax, ai, a2) 

where ps (a^, a^+i, 0-^+2) = Zl^^+a P4 (a^, a^+i, ai+2, a^+s). Full determination of the joint 
probability P4 (aj, a^+i, ai+2, a^+s) requires the calculation of 16 unknowns, but the 
normalization condition and the parity symmetry allows us to reduce this number to 
the following 10 unknowns 

w_i=P4(0,0,0,0), 

Xi=P4(l,0,0,0), X2=P4(0,1,0,0), 

yi=P4(l,l,0,0), 1/2 =P4 (1,0, 1,0), l/s=P4(l,0,0,l), 2/4=^4(0,1,1,0), (43) 

2:1=^4(1,1,1,0), 2:2=^4(1,1,0,1), 

wi =^4(1,1,1,1). 

As usual, Eqs. ([3])- ([6]) allow us to derive the equations for all these unknowns. For 

example, zi = {cricri+iai+2) - (aiai+iai+2ai+s), Wi = (ajai+iai+2ai+s) and so on. 

3.4-1. Mean-field equation Evaluation of the averages in Eqs. ([3])-(|6]) using the 



decomposition (42) results in the following set of equations 
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2 V X2 + yi 



Xi 



4 

2/2 



X2 



yi 

?/3 = - 

1 

^1 = 7 



X2+yi 
X2 



(ys - Xi) + y2 



X2 



yl 



4 \x2 + y2 y2 + Z2J ' 

1 X2y3 

2x2 + yi 



2/2 -- 

?/4 = - 



X2y2 
4 \y2 + Z2 X2+y2 

2/2/2/2 , 2/2 



+ 



4 \X2 + 2/2 2/2 + 2:2 

1 2/4^2 

2 2/1 + ^2 ' 



-a;2 

+ 2 



(44) 



Z2 

2/1 + 2:2 



(2/4 - zi) + 1/2 



1 



4 V3;2 + 2/2 2/2 + Z2 



Wl = — 1 + 

2 V 2/1 + -22 

As the variables w^i and Wi do not affect the other variables and the consistency 
conditions 



'xi ,0-^2 0"a;i ,0'X2 

Y P-ii^,0,a^^,a^^) = Y P4(o-xi, 1,0,(7; 



X2J ) 



X2J 5 



(45 a) 
(456) 



result in the simple relations, 

yi + Z2 = y4 + zi, xi+y3 = X2 + yi, (46) 

from where we can eliminate X2 and yi , we are actually left with a system of 6 coupled 



(47) 



equations for the variables xi, 1/1,7/3, 7/4, 2:1 and Z2. (We note that Eqs. (45a) and (456) 
merely exhibit alternative ways of expressing p2 (1, 1) and p2 (1, 0), respectively.) 
Introducing the linear transformation 

a±=y4±zi, ri± = y2±Z2, 6± = y3±xi 

we obtain the closed set of equations 
1 . 1 

4 . 

r/+ + 4r/+r/_ + r/^ {r]+ + r]_f 



1 



a_ 



V+ + — {V+-V-) 



r/+ = --(27/+ + r/_), ?7_ 

^+ = ^ (a+ + ^- - 5+) , ^- 



m 



a+-V+~ ^-t 



(4J 



5_ 



a+-V+- h + ^ (2a+ - 25+ - r/+ + ?/-; 



Although this system might look formidable, its solution is not very involved. We begin 
by eliminating ?]_ in the equations for a+ and 5+ in order to get 



^ 1. • 1 
c'+ + ^c)+ = a+ + !«+• 

The auxiliary variable / = a+ — 5+ satisfies / = —\f, whose solution is 

/ (i) = («+o - S+o) e--i = po{l- po) (2po - 1) e'i 



(49) 
(50) 
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This explicit solution for / in terms of po and t allows us to consider the equations for 7]± 
as a closed subset of equations which can be solved as follows. The change of variables 



u 



V+ + V- 



7 



V+ + V- 



Vs^ 



v+ f -v+ ' 

leads to the much simpler equations 

1 1 2 1 

uj = --UJ+ -u + g7w, 7'- 

which imply that di^ju) /dt = — |7W. Hence 
7(t)u;(t) = -4po(l-po)e-^ 
where we have used 7 (t = 0) = — 2po and u (t = 0) 



-'JOJ 



(51) 



(52) 



(53) 

2 (1 — Po). Inserting this 
expression back into the equation for 7 we obtain a Ricatti equation, whose exact 
solution is [26] 

7 (t) = -2v/po(l-po)e~3 tanh [S (po, t)] (54) 

with 

S(po,t) = tanh-M^^j +Vpo(l-po)(e"^-l). (55) 

Since we have found explicit solutions for 7 and u we can easily revert to the original 

variables 77+ = f'y/ (7 + u) and //_ = /7 (w — 1) / (7 + u) so as to write 

V+ = -Po (1 - Po) (2po - 1) e-i sinh^ [E (po, t)] (56) 

and 

r?- = -V+ - [po (1 - Po)f^ (2po - 1) e-i sinh [2S (po, t)]. (57) 

At this point we can immediately obtain a_|_, given in Eq. ( [48| ), through a brute-force 
integration 

"+ (^) = oPo (1 - Po) (2po - 1) e-4 



+ 2Po(l 



PoJ 



2 — e 4 1 cosh 



2v/po(l-po) l-e-3 



Po (1 - Po) 



l-2po(l-po)e 3 



sinh 



2VPo(l-Po)(l-e-^)].(58) 



Note that «+ (t — )■ 00) = y (po) calculated within the 3-site approximation [see Eq. 
(30)]. Together with the expression for / given in Eq. (50), this result allows to obtain 
(5+ = «-(. — /. Hence, to complete the solution of the system (48) we need now to 
determine a_ and (5_. This is achieved as follows. 
The auxiliary variable x = ^^ satisfies 

1 



X 



■t(1 + x)^. 
4 a+ 



whose solution is simply 

X(t)=2(l 



Po)e 



-Git) 



(59) 
(60) 
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where we have used x (^ = 0) = 1 — 2po and 

G(i) = j 
Hence a_ (t) = «+ (t) x (t) where the factors in the product are given by Eqs. ( 58 ) 



V+ it') 
it'] 



dt'. 



a_ 



(61) 



and (60). To solve for the last unknown 6- we resort to a shortcut. Considering the 
definitions of 6^ and a;_ in terms of the elementary probabilities given in (43) we see 
that they are related by the transformation 1 -H- and so 6- {t,po) = «- [t, 1 — Po)- 



This concludes the solution of the system ( 48 ) , but we note that since we are not able 



to solve analytically the integral in Eq. (61) the situation here is not as satisfying as 



for the 3-site approximation. Most fortunately, that integral does not appear in the 
expressions for the expectations involving less than 4 sites, as we will see next. 

3.4-2. Calculation of p, (j), ip and Wi. Knowledge of these expectations will allow us 
to compare the 4-site approximation predictions with those of the 3-site approximation. 
The simplest and most important of these expectations is p which can readily be written 
in terms of the previously introduced variables 

p (t) = «+ (t) + 5+ (t) + 7]+ (t) + zi (t) + wi (t) . (62) 

To proceed further we need to derive explicit expression for wi (t). There is a simple 
and elegant way to do that other than replacing the r.h.s. of the equation for wi with 
known quantities and then integrating over t. In fact, introducing e ^ Zi + Wi we have 

^ = ^(2r?+-r7_) (63a) 

= — ?7+ — 2d+, (636) 

(a+ - a_) /2, Z2 = 
yi + Z2 into the equations of ii and wi in system 
(44). As for Eq. {Q3b), it follows directly from the equation for fj+ and d+ in system 
on of e yields 

£(t = 0)-r/+(t = 0) 



where Eq. (63a) was obtained by the direct substitutions Zi 
(r/+-r7_)/2, a^ 



y4 — zi and a+ 



(48). Direct integration of e yields 

e{t) 



2a+ (t = 0) - r/+ it) 
which leads to 

p (t) = £ (t = 0) - 7]+ (t = 0) - 2a+ (t = 0) - / (t) 
= p2 (3 - 2po) + Po (1 - Po) (1 - 2po) e-2 
where we have used / = a+ — 
£ (t = 0) - r/+ (t 



2a+ it) 



5+ given in Eq. (50) and 
= 0) - 2a+ (t 



0)=pa3-2po)=p(po). 



(64) 

(65 a) 
(656) 

(66) 



Most remarkably, Eq. (|656) is identical to its counterpart for the 3-site approximation 



Eq. (23) except for the argument of the exponential in which t/3 is replaced by t/4. 

To derive the 2-site expectation = (crj(Tj+i) we use the relation (j) = yi + Z2 + 
zi + wi = a+ + e so that can be immediately derived using the equations for q;+ 
and 7^+, Eqs. (48) and (56). In this case, the form of the dependence of on po and 
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t has no resemblance with the 3-site approximation counterpart, but the asymptotic 
result is exactly the same. This can be seen by noting that rj^ {t — t- oo) = and so 
= p (po) — a+ (t — !■ oo) which is identical to Eq. (34) since y (po) = a+ (t ^ oo). 

The calculation of the 3-site expectation ip = (crjCri+i(Tj+2) is equally simple. We 
use the relation ip = Zi + Wi = e. Hence ip = p (po) — 2a+ (t — )• oo) which is identical to 



Eq. (35) 



The coincidence between the predictions of the 3-site and 4-site approximations 
for expectations involving 3 contiguous sites provides strong evidence that those 
expectations are the exact results. However, this agreement fails when considering 
expectations involving 4 or more contiguous sites as we can appreciate by calculating 

wi = ((Ti(Ti+icri+20"j+3)- We have wi = e - zi = e - a+ 
and taking t — )■ oo yields 

-G(po) 



[1 - x) /2. Using Eq. (60) for x 



Wiipo) = pipo) + yipo) (1-Po) 



(67) 



where G (po) = f^ [q^ (t') /a+ (t')]dt' and y (po) is given by Eq. (30). As already 
pointed out we have to resort to a numerical integration to evaluate G. Figure [3] shows 
the four expectations calculated in this section. In order to highlight the failure of the 
3-site approximation to estimate 4-site expectations we present in Fig. |4] the comparison 
between the 4-site and the 3-site estimates oiwi. The tiny difference is imperceptible in 
the scale of Fig. [3] but it is sufficient to discard the possibility that 3-site approximation 
is the exact solution of the majority-vote model. 




J 0.0010 

^^ 0.0005 

0.0000 



-0.0005 
0.0 



Figure 4. Discrepancy between tlie steady-state expectation wi ~ {(JiCTi+icr 1+2^^1+5) 



„(4) 



calculated with the 4-site {w\ ) and the 3-site (wj ) approximations 



(3)^ 



3.4-3. Probability of clusters of length m. The procedure here is identical to that 
applied for the 3-site approximation, namely, decompose Pm+2 (1 — o", a, . . . , a, 1 — cr) 
with 0" = 0, 1 in term of the elementary 4-site probabilities. The new element is that 
clusters of length m = 2 can now be described directly by these elementary probabilities. 
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P^PiPo,2) 



y-i + Vi [see Eq. (43)] and yield 

-G(l-po) 



>(4) 



Pr^(po,2)=y(po) 



Poe 



+ (l-po)e 



-GCpo) 



The probability of clusters of length m > 2 is given by 

1 2 / _ \ m— 3 



>(4) 



i^r(Po,m) 






m— 3 



(68) 



(69) 



where ip^i (po) = "^ (1 — Po) and w^i (po) = "U^i (1 — Po) with tp and Wi given by Eqs. 



(35) and (p7j), respectively. These probability distributions are exhibited in Figs, p^ and 
[6j We find a perfect fitting of the simulation data for m = 2; ioi m > 2 the fitting is 
good but there are discrepancies in the vicinity of po = 0.5, which are not perceptible 
in the scale of the figures. 

3.4-4- Two-site correlations- Since the results of the 4-site approximation reduces to 
those of the 3-site for expectations involving up to three contiguous sites, correlations 
such as corr (ctj, ctj+i) and corr (cr.j, crj+2) are the same as for the 3-site approximation. 
In addition {ciCTi+s) ,^) = Wi -\- y^ and so 



(or,or,+3)(4) = P (Po) + V (po) (1 " Po) e'^^"") + poe-'^^^"^'') 



(70) 



However, we need the decomposition in terms of the elementary 4-site probabilities to 
calculate expectations involving more distant sites, such as 



{(^iCi+i) (4) 



W 



i^i;^ 



w_ 



^y 



v^. 



+ 2poy{po)e-^^'-'''\ 



(71) 



These correlations are shown in Figs. [7} Following the already observed pattern, we find 
a perfect agreement with the simulation data for quantities whose calculation involves 
expectations of up to 4 contiguous sites. 

4. Discussion 



Although the main purpose of this contribution is to show the remarkably good 
predictions of the 3 and 4-site approximations to describe the steady-state properties of 
the extended one-dimensional majority-vote model, here we focus on the discussion of 
those properties rather than on the procedure to derive them. 

Figure [5] presents the probability of an absorbing configuration exhibiting a cluster 
of length m as function of the fraction of Is in the random initial configuration. The 3- 
site approximation does not provide a good quantitative account of the simulation data 
but it does provide an excellent qualitative picture which captures the change of Pd 
from unimodal to bimodal that takes place for m = 13. As for the 4-site approximation, 
it provides a very good quantitative representation of the simulation data. In fact, the 
fitting is perfect for m = 2 only, but the discrepancies are so small for m > 2 that 
they are barely visible in the scale of the figure. We note that the transition of the 
distribution Pd {po,m) from unimodal to bimodal was expected. In fact, long clusters 
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Po 



0.00000 Lc 

0.0 0.2 0.4 0.6 

Po 



Figure 5. (Color online) Probability of finding clusters of length m = 2, 13, 15 and 50 
as indicated in the figures. The dashed curves are the results of the 3-site approximation 
and the solid curves of the 4-site approximation. The filled circles are the simulation 
data for a ring of size L — W^ and 10^ independent samples; the error bars are smaller 
than the symbol sizes. The transition from a unimodal to a bimodal distribution takes 
place at m = 13. 



should be abundant for initial configurations with po close to 1 or and very rare when 
the number of Is and Os is well balanced as for po = 0.5. What is surprising is that the 
transition occurs for relatively large m, indicating, for example, that it is more likely to 
find clusters of length 10 by starting with a balanced initial configuration than with an 
unbalanced one. Our simulations indicate that the transition from unimodal to bimodal 
takes place at m = 13 already for L > 17. In addition, the distribution of cluster lengths 
is unimodal for L < 9 for all 2 < m < L — 2. 

Figure |6] shows the dependence of Pd (po, m) on m for fixed po- It is not possible to 
distinguish from the results of the 3 and 4-site approximations and the simulation data 
for Po = 0.2 but some observable discrepancies appear between the 3-site approximation 
and the data for large m in the case po = 0.5. The distribution is given by the sum of 



two exponentials [see Eqs. (38) and (69)] that account for the different possibilities of 
occurrence associated to clusters composed of Is and Os for po ^ 0.5. For po = 0.5 the 
arguments of the two exponentials become identical and so we have a single exponential 
decay. Clearly, for po = 0.2 clusters composed of Is are dominant for small m whereas 
clusters of Os dominate in the large m regime. The slopes of the exponentials are 



complicated functions of po, which can be well-approximated by Eq. (69) derived within 
the 4-site approximation scheme. For comparison. Fig. |6] shows P^i for randomly 
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1x10 



Figure 6. (Color online) Probability of finding clusters of length m for fixed po — 0.2 
(left panel) and po = 0.5 (right panel). The 3-site approximation (dashed curves) 
and 4-site approximations (solid lines connecting the symbols) give results which are 
indistinguishable from the simulation data (filled circles) in the scale of the figure for 
Pq — 0.2; only for very large clusters one can see a noticeable discrepancy between 
the data and the results of the 3-site approximation for pg — 0.5. For the purpose of 
comparison, the solid curves exhibit the results for randomly assembled configurations. 
The simulations were carried out for a ring of size L = 10^ and 10^ independent 
samples. 



assembled configurations which is given by 



-jrandom 
d 



(po,m) = (l-po)V^ + Po(l-Por 



(72) 
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Figure 7. (Color online) Comparison of the results for the 2-site correlations obtained 
with the 3-site approximation (dashed curves), the 4-site approximation (solid curve) 
and the simulations (filled circles). The error bars are smaller than the symbol sizes. 
The left panel exhibits corr {<Ti, (Ti+a) for which the 4-site approximation matches the 
data perfectly, and the right panel shows corr (at, (7^+4) where we can see the failure of 
that approximation scheme. The simulations were carried out for a ring of size L — f 0^ 
and 10^ independent samples. 



The failure of the 4-site approximation in describing all the steady-state properties 
of the extended majority-vote model is better appreciated when we consider the 2-site 
correlations, as shown in FigJ7| As already pointed out, the correlations corr (cij, ctj), 
corr (cri,crj+i) and corr (crj,(Ti+2) are described perfectly by both the 3 and 4-site 
approximations since they involve expectations of two and three contiguous sites only. 
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so Fig. [7] exhibits the more challenging correlations, corr (o"j, ctj+s) (left panel) and 
corr (crj,crj_(.4) (right panel). The 4-site approximation describes perfectly the former 
correlation but not the latter, whereas the 3-site approximation fails in both cases. It is 
interesting that in all cases the 2-site correlations exhibit a peak at po = 0.5. This can 
be explained by noting that the dynamics takes longer to freeze into one of the absorbing 
configurations for the well-balanced initial conditions which results in highly correlated 
sites. On the other hand, for pQ close to its extreme values, most sites are already part 
of frozen random clusters formed during the assemblage of the initial configuration and 
so most of the sites at the final configuration are uncorrelated. 

Figure 7 shows, in addition, that the quality of the approximation improves with 
increasing n, as expected. For example, estimation of corr {ai, crj+4) using the 3-site 
approximation (dashed curve) yields disastrous results, but the results obtained with 
the 4-site approximation (solid curve) are reasonable. We note that as n increases 
the estimation of statistical measures involving n + 1 sites is improved. In fact, the 
relative error resulting from use of the 3-site approximation to estimate corr {ai, ctj+s) at 
Po = 0.5 is 22.5% whereas the error due to use of the 4-site approximation to estimate 
corr (crj,crj+4) is 11.5% only. 

5. Conclusion 

The extended one-dimensional majority-vote model is perhaps the simplest lattice model 
to exhibit an infinity of absorbing configurations. This strong ergodicity breaking is 
probably the reason that the model is not exactly solvable 123 ES]- From the mean-field 
theory perspective, which was the main focus of our paper, the nontrivial nature of the 
steady state of the model presented a most stimulating challenge as the usual fixed-point 
equations proved quite uninformative. In fact, the solution to the problem is a one-to- 
one mapping between the randomly assembled initial configurations, which are described 
statistically by the density po of sites in state 1, and the absorbing configurations. That 
mapping was obtained directly in the case of the pair approximation but in the case of 
the 3 and 4-site approximations we had to solve analytically the dynamics for arbitrary 
t and then take the asymptotic limit t — )■ cxd in order to extract the mapping between 
the initial conditions and the steady state. 

Although the pair approximation describes qualitatively the mapping between po 
and the statistical properties of the steady state, its predictions regarding expectations 
involving two or more contiguous sites are not corroborated by the simulation results 
(see Fig. [2]). The 3-site approximation, however, produces a remarkable good fitting of 
the simulation data for all quantities involving the expectation of three contiguous sites. 
Moreover, the predictions of the 4-site approximation reduce to those of the 3-site in 
the case of three contiguous sites expectations. We see this as a strong indication that 



the expectations {ai), {aiai+2) and (crjCrj+20"j+3) given by Eqs. (24), (34) and (35) are 



exact results. In addition, the perfect fitting of the simulation data by the expectation 



(o"jO"i+20"i+30"i+4)5 calculated within the 4-site approximation [see Eq. (67)] indicates that 
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this quantity may be exact as well, but this evidence is not so strong as for the 3-site 
expectations. 

The findings summarized above as well as a purely numerical analysis of the 5 and 
6-site approximations (data not shown) reveal a most remarkable pattern: the n-site 
approximation seems to yield the exact results for steady-state expectations involving 
n contiguous sites for n > 2. We hope our paper will motivate further research to prove 
(or disprove) this assertion. 
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